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We investigate the suppression of the superconducting transition temperature due to Coulomb 
repulsion in granular metallic systems at large tunneling conductance between the grains, gr ^ 
We find the correction to the superconducting transition temperature for 3-D granular samples and 
films. We demonstrate that depending on the parameters of superconducting grains, the correspond- 
ing granular samples can be divided into two groups: (i) the granular samples that belong to the 
first group may have only insulating or superconducting states at zero temperature depending on 
the bare intergranular tunneling conductance (?t, while (ii) the granular samples that belong to the 
second group in addition have an intermediate metallic phase where superconductivity is suppressed 
while the effects of the Coulomb blockade are not yet strong. 

PACS numbers: 74.81. Bd, 74.78. Na, 73.40.Gk 



A great deal of current experimental and theoretical 
research in mesoscopic physics focuses on properties of 
inhomogeneous granular superconductors 
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interest is motivated not only by their unusual proper- 
ties such as negative magnetoresistance 0, but, even 
to a higher extent, by the fact that they represent an 
exemplary disordered electronic system with experimen- 
tally accessible tunable parameters. A specific feature 
of granular metallic materials is the important role of 
the Coulomb interaction that strongly affects their be- 
havior. The interplay between the Coulomb interaction 
effects and disorder was shown to suppress an electron 
transport in granular metals 0, 0; 13 ; the fundamental 
question that remains open is how the Coulomb repul- 
sion affects the superconducting properties of granular 
metals, and in particular, the superconducting transition 
temperature. 

Coulomb interactions decrease the critical tempera- 
ture, this effect is amplified significantly in the pres- 
ence of disorder and becomes especially pronounced in 
lower dimensions. For example in 2D (disordered) su- 
perconducting films the effect of Coulomb repulsions is 
known to noticeably decrease the critical temperature 
0, Hfl ini H^ . The mechanism of suppression of the 
transition temperature can be understood as a result 
of the renormalization of the electron interaction in the 
Cooper channel by the Coulomb repulsion in the pres- 
ence of scattering by impurities 13]. Since in granular 
samples both disorder and the strong Coulomb repulsion 
effects present, one expects that similar mechanism gov- 
erns the reduction of the superconducting critical tem- 
perate and becomes even more pronounced for granular 
superconductors. 

In this Letter we investigate the influence of Coulomb 
repulsion on the superconducting transition temperature 
of granular metals and demonstrate that depending on 
the parameters of the grains, the granular materials can 
be divided into two groups: (i) materials that may be 
found in either insulating or superconducting state at 



zero temperature and (ii) materials that can, in addition, 
exhibit an intermediate metallic state. The parameter 
that governs the corresponding behavior is the tunneling 
conductance gx- The transport properties of granular 
metals were recently extensively studied in Ref. ]y,|3,|3|. 
It was shown that for 3D samples there is a critical value 
of the bare tunneling conductance 

5? = (1/6^) HEc/S), (1) 

where Ec is the Coulomb energy and 6 is the mean level 
spacing in a single grain, that allows to classify granular 
samples by their transport properties: samples with gx < 
g^ are insulators at zero temperature while samples with 
9t > gT metals. We show that in order to classify 3D 
superconducting granular samples with respect to their 
possible ground state at zero temperature, an additional 
characteristic conductance g^ should be introduced: 

gT^-^n{Ec/T^) (2) 

TT 

where is the superconducting critical temperature of 
a single grain and A is the numerical coefficient defined 
below. We demonstrate that granular superconductors 
can be conveniently classified depending on the relation 
between and g^ into two groups: 

(i) Group A: Granular samples with g^ < gji may have 
only two phases at zero temperature depending on the 
tunneling conductance gx- They are either superconduc- 
tors if gx > gT^ insulators if gx < gx- 

(ii) Group B: Samples with 5^ > g^ can be found in 
three phases. They are superconductors if gx > gxi iii" 
sulators if gx < g^i o^' metals in the intermediate region 

9t> 9t> 9t- 

Thus, in the latter case there exists an intermediate 
metallic region g^ > 9t > 5t where the superconductiv- 
ity is suppressed by the Coulomb interaction while effects 
of Coulomb blockade are not yet strong. 

Presenting our final result for the suppression of the 
superconducting transition temperature we distinguish 



the high temperature, T° > grS, and low temperature, 
T° < qtS, regimes: For 3D systems the suppression of 
the transition temperature is given by 
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while for granular films we obtain 
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is the dimensionless con- 
cos qa) with {a} being 



Here A = gro," 
stant where Eq 
the lattice vectors. Equations Q hold at temperatures 

Turning to the discussion of the results lO, we note 
that the existence of two qualitatively different temper- 
ature regimes (T > qt^ and T < qtS ) is not surprising 
in view of the results for corrections to conductivity ob- 
tained in Ref. where it was shown that the tempera- 
ture dependence of the conductivity at T < is domi- 
nated by the contribution from coherent electron motion 
at large distances, while at T > grS the conductivity 
behavior is controlled by the scales of the order of the 
grain size. From Eq. (|3a|l one can see that the logarith- 
mic corrections appear even in 3D case; this property can 
serve as a specific characteristic of granular metals. In 
the high temperature regime T > similar logarithmic 
corrections (with the same argument but different pow- 
ers) were found in the corrections to the conductivity and 
density of states ^^^^ result for the critical 

temperature suppression in the high temperature regime 
depends on the dimensionality and is not universal con- 
trary to the logarithmic corrections to conductivity found 
in the same regime in Ref. 0, ■ We see that in granular 
superconducting films the suppression of the transition 
temperature is much stronger than that in 31? granular 
samples. The extra logarithm power in two dimensions 
is due to the contribution of low momenta g <C a~^. 

Our classification for 3D granular superconductors can 
be derived from the following arguments: The critical 
conductance is defined by the condition that the cor- 
rection l|Sa|) becomes of the order of unity at gT ~ g^- 
Considering first the case where g^ < gj., we see that 
the system should be a superconductor at gx > g^i 
since the conductance renormalization due to the effect 
of Coulomb interaction l^gx taken at the temperature 
T is smaller than gT for all gT larger than g^. 

Indeed AgT = (I/Gtt) ln(gT-Ec/r°) < c/J,. On the other 
hand, if gT < gT: then the interaction in the Cooper 
channel is completely suppressed and the system becomes 
an insulator at zero temperature because gT < gT- Thus 
in the case 5^ < depending on the value of conduc- 
tance gT one observes either insulating, at gT < gT: or 
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FIG. 1: Diagrams describing the first order self-energy cor- 
rections, 5X1 to the Cooperon channel coupling constant A of 
granular metals due to the electron-electron repulsion. The 
solid lines denote the propagator of electrons, the wave lines 
describe screened Coulomb interaction and the dashed lines 
describe the elastic interaction of electrons with impurities. 
The black rectangle and triangle denote the renormalized 
Cooperon, see Eq. ||KJ, and impurity vertex of granular met- 
als respectively. The indices i and j stands for the number of 
the grains. The tunneling vertices are described by the cir- 
cles. All diagrams should be taken with a factor of 2. The 
final result for the sum of the diagrams (a) - (f) is given by 
Eq. 0. 



superconducting, at gT > gTi phases at zero tempera- 
ture. The second case, where gT > gT: is qualitatively 
different due to appearance of the intermediate metallic 
phase: In the region gT > gT system is a super- 
conductor, in the region gT > 9t > Qt '^'^ interaction 
in the Cooper channel is completely suppressed, and the 
system is metallic since gT > Qt : ^"^^ finally in the region 
9t < 9t fti6 system becomes an insulator. 

Now we turn to the quantitative description of our 
model and derivation of Eq. H3al) : We consider a 
d— dimensional array of superconducting grains in the 
metallic state. The motion of electrons inside the grains 
is diffusive and they can tunnel between grains. We as- 
sume that in the absence of the Coulomb interaction, the 
sample would have been a good metal at T > Tc. 

The system of weakly coupled superconducting grains 
is described by the Hamiltonian 



H — Hq + He + Ht 



(4a) 



The term Hq in Eq. (|4a() describes isolated disordered 
grains with an electron-phonon interaction 

i,k i,k,k' 

(4b) 

where i stands for the number of the grains, k — (k, |), 
— fc = (— k, I); A > is the interaction constant; a|j,(ai,fc) 
are the creation (annihilation) operators for an electron 
in the state k of the i-th grain and Himp describes elastic 
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interaction of electrons with impurities. The term He 
in Eq. 14a() describes the Coulomb repulsion inside and 
between the grains and is given by 



g2 



(4c) 



where Cij is the capacitance matrix and fii is the operator 
of electrons number in the i-th grain. The last term in the 
right hand side of Eq. H4a|l is the tunneling Hamiltonian 



Ht = 



(4d) 



where tij is the tunneling matrix element corresponding 
to the points of contact of i-th and j— th grains. 

The transition temperature, Tc of the granular met- 
als is determined from the pole of the superconducting 
propagator 0| . As in the case with homogeneous super- 
conducting metals in the absence of the electron-electron 
repulsion the transition temperature of granular metals 
is independent of disorder. 

To study the suppression of the transition tempera- 
ture in granular metals we take into account the electron- 
electron repulsion in the Cooper channel. The first or- 
der correction SX to the Cooperon channel coupling con- 
stant A due to electron-electron repulsion is given by 
the diagrams in Figs. 1 and 2. These diagrams con- 
tain essentially the averaged one-particle Green func- 
tions (solid lines), effective screened Coulomb propaga- 
tors (wave lines) and dashed lines that describe the elastic 
interaction of electrons with impurities. In the regime un- 
der consideration all characteristic energies are less than 
Thouless energy Et = D/a? , where D is the diffusion co- 
efficient of a single grain. This allows us to use the zero 
dimensional approximation for a single grain diffusion 
and Cooperon propagators. The electron hopping be- 
tween the grains can be included using the diagrammatic 
technique developed in Refs. 0, 0. The complete expres- 
sion for the renormalized Cooperon of granular metals 
(black rectangle in Fig. 1 and 2) has the following form 



C(w„,q) = (2W) ^{\Lo„\+e^6)- 



(5) 



where q is the quasi-momentum, ujn = 27rTn is the 
bosonic Matsubara frequency and v is the density of 
states on the Fermi surface. The parameter Eq in the 
right hand side of Eq. |(SJl appears due to the electron tun- 
neling from grain to grain, it was defined below Eq. 

Deriving the analytical result for the diagrams in 
Figs. 1 and 2 it is important to take into account the 
fact that the single electron propagator itself gets renor- 
malized due to electron hopping. Tunneling processes 
give rise an additional term to the self-energy part of the 
single electron propagator, see Fig. 3 
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FIG. 2: Diagrams describing the first order vertex corrections, 
5\2 to the Cooperon channel coupling constant A of granular 
metals due to the electron-electron repulsion. AH notations 
are the same as in Fig. 1. The diagrams (b) - (d) should be 
taken with additional factor of 2. After that the final result 
for the sum of the diagrams (a) - (d) is given by Eq. 11111 . 



where tq is the unrenormalized electron mean free time 
and gT = 27rt^i/^ is the dimensionless tunneling conduc- 
tance. Although the second term in the right hand side 
of Eq. © is much smaller than the first one it is impor- 
tant to keep it because the leading order contribution in 
Tg"^ to the coupling constant A cancels ( see Eqs. O and 
^ below). 

We consider the contributions from diagrams in Fig. 1 
and 2 separately. Using the result of Eq. ^ it is straight- 
forward to check that the sum of the diagrams (a)-(f) 
in Fig. 1 results in the following correction, 5\i to the 
Cooperon channel coupling constant 
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Here the summation is going over the quasi-momentum, 
q and over the fermionic, e„ — 7rT{2n + 1) and bosonic, 
fin — 2'iTTn Matsubara frequencies. The propagator of 
the screened electron-electron interaction, y(ri„,q) in 
Eq. Q is given by the expression 



^(f^„,q) = 



2gc(q)(|»„|-t-£qJ) 
4eq£:c(q) |f^„| 



(8) 



The parameter gq in Eqs. (TJ and ijHJl was defined below 
Eq. H3a|) . The charging energy Ec{q) = e'^/2C{q) in 
Eq. ((SJ is expressed in terms of the Fourier transform 
of the capacitance matrix C(q) which has the following 
asymptotic form 



C-'i^) - 777 



■n/q 



D = 2, 
£> = 3. 



(9) 
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FIG. 3: Diagrams describing self-energy corrections to the 
single electron propagator due to a) elastic interaction of elec- 
trons with impurities and b) electron hopping. The solid lines 
denote the bare propagator of electrons and the dashed line 
describes the elastic interaction of electrons with impurities. 
The tunneling vertices are described by the circles. 



To obtain the total correction, SX to the Cooperon 
channel coupling constant A due to electron-electron in- 
teraction the diagrams in Fig. 2 that represent vertex cor- 
rections should be taken into account. These diagrams 
results in the following contribution 



is suppressed. Taking into account effects of Coulomb 
blockade of conductivity that become essential for sam- 
ples with 5t < 5t ! where the critical value is given by 
Eq. we introduced the classification of 3-D granular 
samples composed from superconducting grains accord- 
ing to the possibility of having insulating, metallic or 
superconducting phases at zero temperature. 
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Using Eqs. ||SJ|, lO and summing over the frequencies e„, 
fin and quasi-momentum, q in Eqs. |7J and H10|) with 
the logarithmic accuracy we obtain the final result (jJJ for 
the superconducting transition temperature in granular 
metals. 

In conclusion, we considered the suppression of super- 
conducting transition temperature due to Coulomb re- 
pulsion in granular metals. We found the correction to 
the transition temperature for 3-D granular samples and 
films at large tunneling conductance between the grains. 
We demonstrated that the suppression of superconduc- 
tivity in 3-D samples can be characterized by the critical 
value of the conductance introduced in Eq. ^ such 
that for samples with gx < gx the superconductivity 
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